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ON  THE  CONSERVATION  AND  CONVERGENCE  TO  WEAK  SOLUTIONS  OF 

GLOBAL  SCHEMES* 

MARK  H.  CARPENTERt,  DAVID  GOTTLIEBt,  AND  CHI- WANG  SHU§ 


Abstract.  In  this  paper  we  discuss  the  issue  of  conservation  and  convergence  to  weak  solutions  of  several 
global  schemes,  including  the  commonly  used  compact  schemes  and  spectral  collocation  schemes,  for  solving 
hyperbolic  conservation  laws.  It  is  shown  that  such  schemes,  if  convergent  boundedly  almost  everywhere, 
will  converge  to  weak  solutions.  The  results  are  extensions  of  the  classical  Lax-Wendroff  theorem  concerning 
conservative  schemes. 

Key  words,  conservation  laws,  conservation,  weak  solutions,  convergence 
Subject  classification.  Applied  and  Numerical  Mathematics 

1.  Introduction.  We  are  interested  in  numerical  solutions  to  the  conservation  laws: 


(1.1) 


Ut  +  f{u)x  =  0,  M(a:,0)  =  u^{x),  — 1  <  a;  <  1. 


Here  we  have  written  (1.1)  in  the  one  dimensional  form,  but  the  results  of  this  paper  are  also  valid  for  multi 
dimensions. 

The  purpose  of  this  paper  is  not  to  study  the  issue  of  convergence.  We  actually  assume  that  the  numerical 
solution  converges  boundedly  a.e.  (almost  everywhere),  to  a  certain  function  u{x,t).  More  precisely,  for  a 
numerical  scheme  defined  at  the  (uniform  or  nonuniform)  grid  points  Xj,0  <j  <  N,  with  Ax  =  max(a;j+i  — 
Xj)  and  Vj{t)  as  the  numerical  solution  at  x  =  Xj,  we  define  the  function  VAx{x,t)  by 

(1-2)  VAx{x,t)  =  Vj{t),  Xj  <X  <  Xj+i, 


and  assume  that  VAx{x,t)  is  uniformly  bounded  with  respect  to  x,  t,  and  Ax,  and,  as  Ax  -t  0,  VAx(x,t) 
converges  pointwise  a.e.  to  u(x,t).  See,  e.g.,  [5,  16, 17,  20]  for  discussions,  in  the  scalar  case,  of  convergence 
of  some  of  the  schemes  studied  in  this  paper,  under  the  L°°  boundedness  assumption.  We  will  concentrate 
on  the  issue  of  whether  the  limit  function  u(x,t)  is  a  weak  solution  to  (1.1),  that  is  whether  it  satisfies 
rT 

(1.3) 


-L 


t)4>t{x,t)  f{u{x,t))<f>x{x,t))  dxdt—  I  u^{x)(l}{x,0)dx  =  0 


i: 


for  any  smooth  function  <f){x,  t)  which  is  compactly  supported.  Also,  in  this  paper  we  only  consider  semidis¬ 
crete  method-of-lines  schemes,  i.e.  schemes  which  are  discretized  in  the  spatial  variable(s)  only. 

The  classical  result  in  this  area  is  the  famous  Lax-Wendroff  Theorem  [11]: 


Theorem  1.1.  (Lax  and  Wendroff)  If  the  numerical  solution  of  a  conservative  scheme: 
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where  the  numerical  flux 

(1-5)  fj^L=f{Vj-p,...,Vj+q) 

is  local  (i.e.  p  and  q  are  constants  independent  of  Ax),  Lipschitz  continuous  in  every  argument,  and  consistent 
with  the  physical  flux  f{v,  ...,v)  —  f{v),  converges  boundedly  a.e.  (almost  everywhere)  to  a  function  u{x,  t), 
then  u{x,t)  is  a  weak  solution  to  (1.1). 

■ 

The  proof  follows  easily  from  a  summation  by  parts  and  an  application  of  the  dominant  convergence 
theorem.  See  LeVeque  [13]  for  a  slightly  different  version  of  this  theorem  and  its  proof. 

The  Lax-Wendroff  Theorem,  however,  does  not  cover  global  schemes,  i.e.  schemes  which  can  not  be 
written  in  the  form  (1.4)  with  a  local  flux  /j+i .  Excimples  of  global  schemes  include  the  compact  schemes 
[12,  2,  3,  5],  and  various  spectral  Galerkin  or  collocation  schemes  (Fourier,  Legendre,  Chebyshev)  [9,  1,  6]. 
We  will  extend  the  Lax  Wendroff  Theorem  to  these  global  schemes  in  this  paper. 

We  remark  that  there  are  discussions  in  the  literature  about  schemes  which  are  not  of  the  conservative 
form  (1.1)  but  nevertheless  still  converge  to  weak  solutions.  One  such  example  is  the  class  of  schemes  for 
general  curvilinear  coordinates,  see  [19]  for  a  proof  that  such  schemes  actually  do  converge  to  weak  solutions. 

In  [8]  the  authors  discussed  conservation  issues  of  Chebyshev  methods.  However,  they  only  considered 
the  mean  of  the  solution,  that  is,  they  verified  that  the  limit  solution  satisfies  (1.3)  with  <f>{x,  t)  =  1. 

The  organization  of  the  paper  is  as  follows.  In  Section  2  we  discuss  compact  schemes.  Section  3  contains 
the  Legendre  collocation  method,  while  Section  4  discusses  the  Chebyshev  method.  Section  5  discusses  the 
Legendre  approximation  in  the  multi-domain  case. 

To  close  this  section  we  mention  that  in  this  paper  C  (or  c)  is  a  generic  constant. 

2.  Compact  Schemes.  Compact  schemes  are  methods  where  the  derivatives  are  approximated  by 
rational  function  operators  on  the  discrete  solutions.  We  consider  compact  schemes  defined  on  a  uniform 
grid  Xj,0  <j<N.  For  example,  a  fourth  order  central  compact  approximation  to  the  derivative  is  [12]: 

(2.1)  i  ((«x)i-i  +  +  (ux)j+i)  =  ^  {vj+i  -  Vj-i) 

and  a  third  order  upwind  compact  approximation  to  the  derivative  is  [5]: 

(2.2)  -  +  h{Vx)j  —  (nx)j+i)  =  —  (3vj  —  4Uj_i  +  Vj-2)  ■ 

Adequate  boundary  conditions  must  be  used  for  the  compact  schemes,  to  retain  accuracy  and  stability,  see 
[2],  [3]  for  details.  Together  with  boundary  conditions,  a  compact  scheme  for  (1.1)  can  be  written  as 

(2.3)  Pvt  -I-  Qf{v)  -  -t{vb  -  9b) 

where  v  =  (uo,  ■■■Vn)*  is  the  numerical  solution,  t  is  a  constant,  vb  =  (uojO,  ...,0,U7v)*  is  the  boundary  part 
of  the  numerical  solution,  and  qb  =  (go,0,  ...,0,9jv)*  is  the  given  boundary  data.  Depending  on  the  wind 
direction,  one  or  both  of  the  first  and  last  components  of  ub  and  qb  may  also  be  zero(s).  The  matrices  P 
and  Q  satisfy  the  following  conditions  [3]: 

•  P  is  symmetric,  and  satisfies 

(2.4)  P4>-(j)  =  0{Ax) 

Here  and  below,  (f>  =  {(p{xo), ...,  i^(xjv))*,  and  (j){x)  is  an  arbitrary  smooth  (C^  or  smoother)  function. 
0{Ax)  for  a  vector  means  that  each  component  is  bounded  by  a  constant  times  Ax,  and  the  constant 
depends  only  on  the  derivatives  of  <j>{x); 
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Q  is  “almost”  anti-symmetric,  that  is: 


Q  +  Q*  =  R  +  S 


where  R  =  (ry),  and  rij  =  0  except  for  roo  and  r^N.  5  is  either  identically  0  for  the  central  compact 
schemes,  or  satisfies 

(2.6)  S(l>  =  OiAx) 

for  the  upwind  compact  schemes,  where  (p  is  defined  as  before.  Also,  Q  is  at  least  a  first  order 
approximation  to  the  derivative: 

(2.7)  Q4>-<t>^=OiAx). 

We  can  easily  verify  that  all  the  compact  schemes  in  [12,  2,  3,  5]  satisfy  the  above  conditions  for  P  and 

Q- 

For  such  compact  schemes  we  can  state  the  following  proposition: 

Proposition  2.1.  If  the  solution  of  the  compact  scheme  (2.3)  converges  almost  everywhere  to  a  function 
u{x,t),  then  u{x,t)  is  a  weak  solution  to  (1.1). 

Proof:  For  any  compactly  supported,  function  <f>{x,t),  we  denote  by  if>  =  {(j>{xo,t),  ,  left 

multiply  (2.3)  by  (p*,  and  integrate  over  [0,T]  to  obtain: 

f  <p*{Pvt  +  Qf(v))  dt  =  0 
Jo 

due  to  the  zero  boundary  conditions  of  </>.  Now  integrating  by  parts  in  t  for  the  first  term,  taking  a  transpose 
of  the  equation  (which  is  a  scalar),  and  using  the  symmetry  of  P  and  condition  (2.5)  of  Q,  we  obtain: 

-  r  {v*P<Pt  +  f{vyQ<p)dt-{v*P<P)\t=o  =  -  r  f{vyS(Pdt 
Jo  Jo 

Or,  considering  (2.4),  (2.6),  (2.7),  and  the  uniform  boundedness  (with  respect  to  the  mesh  size  Ax)  of  v, 

(2.8)  -  r  {v*<Pt  +  f{vy<py,  dt  -  {v*<p)  l«=o  =  0(1) 

Jo 

where  the  constant  term  0(1)  results  from  a  summation  of  N  =  terms  of  0(Ax)  quantities;  the 

constant  depends  only  on  the  derivatives  of  (p{x). 

Recalling  the  definition  of  the  function  V£^x{x,t)  in  (1.2),  we  can  multiply  (2.8)  by  Ax  to  obtain 

{2.Q)  -  J  j  [vAx{x,t)4>t'^{x,t)  +  f{v^x{x,t))<Px^(x,t))dxdt-J  VAx{xyO)(p^^  {x,0)  dx  =  0{Ax) 


(2.10)  <p^^{x,t)  =  (p{Xj,t),  ^^^{x,t)  =  (px{Xj,t),  <pf'^{x,t)  =  (pt{Xj,t),  Xj  <  X  <  Xj+i, 

By  assumption,  VAx{x,t)  converges  to  u(x,t)  boundedly  a.e.  There  is  no  problem  about  the  uniform 
convergence  of  <p^^(x,t),  0^^(x,t)  and  4>^^(x,t)  due  to  the  smoothness  of  0.  By  the  dominant  convergence 
theorem,  taking  the  limit  as  Ax  — >  0  in  (2.9),  we  obtain  (1.3).  This  proves  that  u(x,t)  is  a  weak  solution  of 
(1.1). 
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3.  Legendre  Spectral  Collocation  Schemes.  The  Legendre  collocation  method  can  be  written  in 
the  following  way: 

where  un{x,  t)  is  the  numerical  solution  which  is  a  polynomial  of  degree  at  most  N  in  x,  In  is  the  Legendre 
interpolation  operator,  i.e.  for  any  function  g{x),  1n9(x)  is  the  unique  polynomial  of  degree  at  most  N 
satisfying  lNg{xj)  =  g{xj)  at  the  +  1  Legendre  Gauss- Lobatto  points  xj,  which  are  the  zeros  of  the 
polynomial  (1  —  x^)Pj^i,  where  Pn  is  the  Legendre  polynomial  of  degree  N. 

The  term  SV  is  the  spectral  viscosity  term  needed  to  stabilize  the  scheme  and  in  order  for  the  assumption 
“vax(x,  t)  converges  boundedly  a.e.  to  a  function  u(x,  t)”  to  be  realistic.  We  consider  here  the  superviscosity 
term 

(3.2) 

£  is  the  superviscosity  coefficient,  s  is  an  integer  growing  with  N  [21,  10,  14,  15].  We  remark  that  this 
superviscosity  term  is  equivalent  in  practice  to  a  low  pass  filter. 

Finally,  the  boundary  term  Bujv(a:,  f)  could  be  either  0,  or 

T(uAr(l,t)  -5i)(l  +  a:)P;,f(x), 

or 

T(uN(-l,t)  -  g-i)(l  -  x)Pn(x), 

or  a  combination,  depending  on  the  wind  directions  at  the  boundary  points.  Here  t  is  a  constant  chosen  for 
stability  and  .91  and  9_i  are  functions  of  the  time  only. 

Let  (f>{x,t)  be  a  test  function  in  C^.  Take<^;v-i(x,  t)  =  lN-\4>{x,t),  then  clearly  0jv-i  are  polynomials  of 
degree  at  most  W-1  and  vanish  at  both  boundary  points  x  =  ±1.  Also  (t>N-i{x,  t)  ->  cf>{x,  t),  {4>n-i)x{x,  t) 
4>x{x,t),  and  {4>N-\)t{x,t)  (j>t{x,t)  uniformly. 

We  denote  now  by 

(/,<?)  =  j  J{x)g{x)dx 

and  by 

N 

{f,9)N  =  Y^f(xj)g{xj)u)j 
k=0 

where  luj  >  0  are  the  weights  in  the  Gauss-Lobatto  formula.  We  note  that  {f,g)N  —  if >9)  if  f9  is  a 
polynomial  of  degree  at  most  2N  —  1. 

We  first  show  that  the  boundary  terms  do  not  cause  a  problem: 

Lemma  3.1. 

(3.3)  jBmjv)  =  0. 
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Proof:  We  start  by  observing  that 

{(I>N-1,Bun)  =  i<j>N-l,BuN)N. 

Bun  vanishes  for  the  inner  Gauss-Lobatto  points  and  (j>N-i  vanishes  at  the  boundaries  and  therefore  the 
Lemma  is  proven. 


With  (3.3),  we  multiply  (3.1)  by  cj)N-i{x,t),  integrate  over  x,  and  integrate  by  parts  for  the  second  term 
to  obtain 

J  (j>N-\ix,t) - — - ~  J  - - lNf{uN{x,t))dx 

(3.4)  =  UN{x,t)dx. 

We  now  estimate  the  right  hand  side  of  (3.4): 

Lemma  3.2. 

(3.5)  jim^  4>N-i{x,t)  [^(1  UN{x,t)dx  =  0. 

Also,  the  quantity  under  the  limit  sign  is  uniformly  bounded  with  respect  to  t. 

Proof:  Since  (f>N-i  is  a  polynomial  of  degree  AT  —  1, 

N-l 

(l)N-l{x,t)  =  V  ^k,N{t)Pkix) 
k=0 

where  ^k,N(t)  are  the  collocation  Legendre  coefficients  of  the  test  function  (j).  Note  that 

^(1  -  x^)^y  Pk{x)  =  i-iyk^k  +  lyPkix), 

and  therefore 


N 

^  k%k  +  iy$k,N(t)uk,N{t){Pk,Pk), 


k=0 


Here  Uk,N  are  the  Legendre  collocation  coefficients  of  un-  We  note  that  as  a  consequence  of  the  uniform 
boundedness  of  UN{xj)  and  the  fact  that  is  in  Cq°, 

(3-6)  \^k,N(t)Uk,N{t)\  < 

This  implies  (3.5)  and  the  uniform  boundedness  of  the  quantity  under  the  limit  sign  with  respect  to  t. 


We  thus  only  have  to  deal  with  the  left  hand  side  of  (3.4).  We  integrate  (3.4)  in  t,  integrate  by  parts  for 
the  first  term,  and  use  Lemma  3.2  to  obtain: 


fT  pi  / 

■  /  /  ( UN{x,i 

Jo  J-i  \ 


j  ur^{x,0}(j>r^-i{x,G)  dx  —  o(l). 
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It  looks  like  we  can  immediately  take  the  limit  as  in  the  Lax-Wendroff  Theorem.  The  trouble  is  that,  we 
have  only  assumed  the  uniform  boundedness  of  Ufj{xj,i),  hence  of  f{uNixj,t)),  but  this  does  not  imply  the 
uniform  boundedness  of  either  UN{x,t)  or  /w/(wjv(a;,  t))  due  to  the  lack  of  regularity. 

We  need  the  following  Lemma; 

Lemma  3.3.  Let  vax  be  the  piecewise  linear  polynomial  taking  the  values  UN{xj,t),  then 

(3.8)  =  +<.{!), 

(3-9)  [iNtMt,  i)),  +0(1). 


Proof:  We  will  switch  back  and  forth  between  integrals  and  quadrature  summations: 


U/\/  ^Xy  t)  , 


d^N--i{x,t) 

dt 


UN{x,t), 


d(j>N-\{x,  t) 
dt 


N 


N 


^\d(j)N-\{Xj,t) 
2_^UN(Xj,t) - - LJi 


i=o 


dt 


=  j  e{x)vAx{x,t){(j>N-l)f'''{x,t)dx 


where 


6(x)  = - — - ,  Xj  <x  <  Xj+i. 

Xj+l  -  Xj 

In  [18]  it  has  been  established  that  d(x)  is  uniformly  bounded  and  converges  a.e  to  1  as  A  a;  0. 
This  proves  (3.8).  The  proof  for  (3.9)  is  similar. 


We  can  state  now 


Theorem  3.4.  If  the  function  VAx{x,t)  defined  in  (1.2),  obtained  from  the  solution  of  the  Legendre  collo¬ 
cation  scheme  (3.1)  at  the  Legendre  Gauss- Lobatto  points  xj,  converges  almost  everywhere  to  a  function 
u{x,t),  then  u{x,t)  is  a  weak  solution  to  (1.1). 

Proof:  By  assumption,  VAx{x,t)  converges  to  u{x,t)  boundedly  a.e.  Also,  there  is  no  problem  about  the 
uniform  convergence  of  (f>%^_-^{x,t),  {4>N-{)x^{x,t)  and  {<pAi-i)^^(x,t)  due  to  the  smoothness  of  (j>.  Using 
the  dominant  convergence  theorem,  taking  the  limit  as  Ax  — »•  0  in  (3.8),  (3.9)  and  (3.7),  we  obtain  (1.3). 
This  proves  that  u{x,t)  is  a  weak  solution  of  (1.1). 


We  close  this  section  by  commenting  on  other  spectral  viscosity  terms  in  (3.1)  that  stabilize  the  Legendre 
method.  One  such  term  is 


cjv 


d  „duN 
dx  dx 


where  the  spectral  viscosity  operator  Q  is  defined  by 


N 

Q<i>  =  ^<3ifc<^ifcP/fc(a:) 

k=0 
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where 


N 


(j> ='Y^^kPk{x) 

it=0 


and 


Qk  =  0,k  <  tun 
1  >  <3*  >  1  -  (^) 

with  mAT  growing  with  N. 

We  can  establish  also  for  this  viscosity  term  that 


,  SV (un))  0 


and  therefore  the  result  above  holds  also  for  this  kind  of  spectral  viscosity. 


4.  Chebyshev  Spectral  Collocation  Schemes.  In  this  section  we  consider  the  Chebyshev  collocation 
schemes.  These  are  more  difficult  to  analyze  than  the  Legendre  method  because  of  the  weight  function  ■ 
The  Chebyshev  collocation  method  can  be  written  in  the  following  way: 


(4.1) 


duN{x,t)  dJNf{UN{x,t)) 
dt  dx 


UAr(x,t)  +  BuN{x,t) 


where  again  UN{x,t)  is  the  numerical  solution  which  is  a  polynomial  of  degree  at  most  N  in  x,  Jn  is  the 
Chebyshev  interpolation  operator,  i.e.  for  any  function  g(x),  Jn9{x)  is  the  unique  polynomial  of  degree  at 
most  N  satisfying  Jffg{xj)  —  g{xj)  at  the  N  +l  Chebyshev  Gaus.s-Lobatto  points  Xj.  e  is  the  superviscosity 
coefficient,  s  is  an  integer  growing  with  N  [21, 14, 15].  We  remark  again  that  this  superviscosity  term,  which 
in  practice  is  equivalent  to  a  low  pass  filter,  or  a  similar  vanishing  viscosity  term  [16,  17],  is  needed  in  order 
for  the  assumption  “vAx{x,t)  converges  boundedly  a.e.  to  a  function  u{x,ty’  to  be  realistic.  Finally,  the 
boundary  term  Bun{x,  t)  could  be  either  0,  or 


t{un{1,  t)  -  5i)(l  +  x)T}^{x), 


or 


r(ujv(-l,t)  -  (/-i)(l  -  x)T'^{x), 

or  a  combination,  depending  on  the  wind  directions  at  the  boundary  points.  Here  r  is  a  constant  chosen  for 
stability  and  gi  and  g-i  are  functions  of  the  time  only. 

Let  0(a;,  t)  be  a  test  function  in  Cq,  that  is,  all  x  derivatives  of  (j){x,  t)  up  to  order  5  vanish  at  the  boundary 
points  X  =  ±1.  Such  test  functions  are,  of  course,  dense  in  Cq.  It  follows  that  (1  -  is  in  Cq. 

We  denote  the  (JV  —  5)-th  degree  Chebyshev  interpolation  polynomial  of  the  function  (1  —  x^)~^^^(j>{x,  t)  by 

(4.2)  ^N-Q{x,t)  =  Jjv-5((1  -  x'^y^l‘^<i>{x,t)), 

and  note  that 


^N-5ix,t)  (1  -x^) 

d^N^5{x,t)  d  {{1  -  x'^)-^/'^(f>{x,t))  O^nMx,  t)  d  {{I  -  x^)-^/‘^(t>{x,t)) 
dx  dx  ’  dt  dt 
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uniformly.  We  now  take 


ll)N-l{x,t)  =  (1  - 


then  j/’jv-i  is  a  polynomial  of  degree  at  most  iV  —  1  and  vanishes  at  both  boundary  points  a;  =  ±1  together 
with  its  first  and  second  x  derivatives.  Moreover,  it  can  be  easily  verified  that 


uniformly. 

We  again  first  show  that  the  boundary  term  does  not  cause  a  problem: 


Lemma  4.1. 


Ijl  +  x)r^{x)'^^^^dx  =  0,  - x)T'„ix) 


>f  /  \  1  (^)  0 


dx  =  0. 


Proof:  We  only  prove  the  first  equality.  Zero  boundary  values  of  ipi\i-i{x,  t)  and  its  first  x  derivative  imply 

J\l  +  x)T^ix)'t^^^  dx  =  -  £  T^(a:)  +  a:))  dx 


C  Ill^( 

7-1  \/l  -  x'^  V 


d’lpN-\{x,t)  ,,  ,  ,  'tjjN-l(x,t) 

— aj — (1  +  “^)+ 


The  last  equality  is  due  to  the  fact  that 

is  a  polynomial  of  degree  at  most  W  —  1,  hence  is  orthogonal  to  Tn{x)  with  the  weight  :^7j=^- 

■ 

With  (4.5),  we  can  now  multiply  (4.1)  by  integrate  over  x,  and  integrate  by  parts  for  the 

second  term  to  obtain 


^  IpN-iiXjt)  duN{x,i) 


r 

7-1 


e2  dt 


!ldx-  f 

7-1  \  %/i  — 

)^  /■*  ^Ar-i(3:,t)  r  _ 

-i7-,  ■ 


JNf(uN{x,t))dx 


a;^)—  UN{x,t)dx. 


We  now  estimate  the  right  hand  side  of  (4.6): 


Lemma  4.2. 


,i„  ^  f 

N^oo  JV2s  1  J_^  y/l  — 


'^TT-  UN{x,t)dx  —  0. 
ox 


Also,  the  quantity  under  the  limit  sign  is  uniformly  bounded  with  respect  to  t. 

Proof:  Integrating  by  parts  2s  times,  and  noticing  that  the  boundary  terms  are  always  0  because  of  the  fact 
that  '4>N-i{x,t)  vanishes  at  the  boundaries  and  because  of  the  factor  \/l  —  a:^,  we  obtain 


<“■*>  L 


’uN(x,t)dx  =  y' 
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Recalling  the  definition  of  ^N-&{Xyt)  in  (4.2),  we  have 


^N-s{x,t)  =  ^  ik{t)Tk{x) 


where  ^k{t)  are  the  collocation  Chebyshev  coefficients  of  the  Cq  function  (1  —  x^)  ^/‘^(f){x,t),  hence 
(4.9)  \Ut)\  < 


Now,  by  the  relationship  between  and  ^at-s  in  (4.3): 


ipN-i{x,t)  =  (1  -x^)^  ^  ik(t)Tk{x) 


=  \{l-Ux)r'^Ut)Tk{x) 

k=0 

-  iV-5  N-5 

=  4  E  Ut)Tk{x)  --J2  Ut)  {Tk+2ix)  +  TkM^)) 

^  k=0  ^  Jfc=0 

,  N-5 

+  T^  E  (^^+''("^)  +  2n(x)  +  Tk-4(x)) 

°  k=0 

“  16  E  -  •iik-2{t) + 6|fc(t)  -  iik+2{t) + ik+i{t)j  Tk{x) 

k=Q 

N~\ 

=  ^  Mt)Tk{x) 

Jfe=0 

where  we  take  the  convention  that  ^^(t)  =  0  for  A;  <  0  or  ^  —  5.  This,  together  with  (4.9),  clearly 

implies 


We  now  use  the  equality 


\M)\  < 


il’N-i{x,t)  =  '^{-iyk‘^^'ipk{t)Tk{x) 


and  the  integral-quadrature  equivalence: 


I  I  1pN-li,X,t)dx  =  ^WjUN{Xj,t)  ^ 


■a;2—  ipN-i{x,t) 


where  Xj  and  Wj  are  the  nodes  and  weights  of  the  Chebyshev  Gauss-Lobatto  quadrature  formula,  because 
the  integrand 


UNix,t)  ^N-iix,t) 

is  a  polynomial  of  degree  at  most  2JV  —  1.  This,  together  with  the  uniform  boundedness  of  ujvixjjt)  and 
(4.8),  (4.10)  and  (4.11),  implies  (4.7)  and  the  uniform  boundedness  of  the  quantity  under  the  limit  sign  with 
respect  to  t. 
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We  thus  only  have  to  deal  with  the  left  hand  side  of  (4.6).  We  integrate  (4.6)  in  t,  integrate  by  parts  for 
the  first  term,  and  use  Lemma  4.2  to  obtain: 


t)) 


(4.12) 


L  I  1  +JNf(uNix, 

-  j  un{x,Q) 


y/\  —  Xp' 

V-'jv-i(z,0) 
%/!  — 


dxdt 


dx  —  o(l). 


Again,  the  difficulty  is  that  we  have  only  assumed  the  uniform  boundedness  of  ui^{xj,t),  hence  of  f{uN{xj,t)), 
not  the  uniform  boundedness  of  either  Uf^{x,i)  or  J]\jf{ui\[{x,t)).  We  again  get  around  this  by  switching 
between  integrals  and  quadrature  summations: 


(4.13) 


(W). 

==  J2wj\/l-X^  (uNiXj,t) 


+  JNf{uN{x,t)) 


'4>N-\{x,t) 

\/l  —  X^ 


dx 


+  f{uN{Xj,t)) 


/  1pN-l{x,t)\ 

V  v'l  -  a;2  ) , 


X  —  Xj 


because  we  can  easily  verify  that  the  integrand 

is  a  polynomial  of  degree  at  most  2N  —  1.  Recalling  the  definition  of  the  function  Vj\x{x,  t)  in  (1.2)  and  that 
of  4>^^{x,t)  etc.  in  (2.10),  we  can  use  (4.13)  to  rewrite  (4.12)  as 


(4.14) 


-  j  0(x)uAi(a:,O) 


t)) 


tpN-i{x,t) 

\/l  —  X^ 


Ax 


dxdt 


^N-l{x,0) 

V^l  -  x^ 


Ax 


dx  =  o(l) 


where 


0{x) 


Wj\/1  —  x'^ 
Xj-j-l  ■  Xj 


Xj  <X  <  Xj^i . 


Clearly,  0{x)  is  uniformly  bounded  and  converges  to  1  as  Aa;  ->  0.  By  assumption,  VAx.{^-i)  converges  to 
u{x,  t)  boundedly  a.e.  Also,  (4.4)  guarantees  the  uniform  convergence  of  the  ipN-i  related  terms  to  the  right 
limits.  Using  the  dominant  convergence  theorem,  taking  the  limit  as  As  ->  0  in  (4.14),  we  obtain  (1.3). 
This  proves  that  u{x,t)  is  a  weak  solution  of  (1.1),  i.e.  we  have  proved  the  following 


Proposition  4.3.  If  the  function  VAx{x,t)  defined  in  (1.2),  obtained  from  the  solution  of  the  Chebyshev 
collocation  scheme  (4.1)  at  the  Chebyshev  Gauss-Lobatto  points  xj,  converges  almost  everywhere  to  a 
function  u{x,t),  then  u{x,t)  is  a  -weak  solution  to  (1.1). 


5.  Multi-Domain  Legendre  Methods.  In  this  section  we  will  discuss  stable  and  conservative  inter¬ 
face  boundary  conditions  for  the  multi-domain  Legendre  method  applied  to  equation  (1.1).  We  assume  that 
the  domain  -1  <  s  <  1  is  divided  into  two  domains,  and  for  the  sake  of  simplicity  we  assume  that  the 
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interface  point  is  a;  =  0.  We  will  denote  by  UN(x,t)  the  numerical  approximation  in  — 1  <  x  <  0  and  by 
VN{x,t)  the  solution  at  0  <  x  <  1.  The  multi-domain  Legendre  method  is  given  by 

~  B{uN{-\,t))  +  nQ‘ {x)  [/+(ujv(0,t))  -/+(ujv(0,t))] 

(5.1)  +T2<3^(x)  [/“(ujv(0,t))  -  f^{vN{Q,t))]  +  SViuN), 

^  nQ^ix)  [/+(uN(o,t))  -/+(u7v(o,t))] 

+TiQ^'{x)  [f~ivN(0,t))  -  f~iuN{0,t))] 

(5.2)  +Sn«Jv)+B(r;jv(l,t)). 


Equation  (5.1)  holds  in  the  interval  — 1  <  a;  <  0,  and  (5.2)  holds  in  0  <  x  <  1.  lj^f{uN)  interpolates  /(uat) 
at  the  zeroes  of  the  polynomial  xQ^  and  Ijj  f{vN)  interpolates  /(ujv)  at  the  zeroes  rjj  of  the  polynomial 
xQ^^ ,  where 


(1  +  x)p;(2x  + 1)  ,,  _  (i-x)p;(2x-i) 

Pa-(1)  ^  ’  ^a(-I) 


The  spectral  viscosities  SV{un)  and  SV{vn)  are  of  the  form 


(5.3) 

(5.4) 


(-11®+’ 


UN, 


VN- 


At  this  point  we  stress  that  the  results  of  this  section  are  valid  only  for  this  form  of  spectral  viscosity 
and  not  for  the  others  discussed  in  Section  3.  The  reason  for  that  will  be  evident  in  the  proof. 

Finally  the  boundary  operators  B  at  the  ends  of  the  interval  — 1  <  x  <  1  are  left  unspecified  for  now. 

We  will  also  denote  the  scalar  product  {p,q)N  =  P(®)  ?(*)  defined  in 

[-1,0]  and  {p,q)N  =  Z]j=oP^(’lj)9(^i)‘^i  f(^)  9(^)  defined  in  [0,1],  and  Wj  are  the  weights  in 

the  Gauss-Lobatto  Legendre  quadrature  formula.  Note  that  if  pq  is  a  polynomial  of  degree  at  most  2N  -  1 
defined  in  [—1,0]  then 

fO 

ip,g)N  =  {p,q)  =  J  ^p^{x)q{x)dx. 


A  similar  formula  holds  in  the  interval  [0, 1]. 

Our  aim  in  this  section  is  to  show  that  the  choice  of  the  parameters  Tj,  i  =  1,4  that  leads  to  linear 
stability  is  sufficient  for  proving  conservation,  i.e.  if  the  numerical  solution  UN{,x,t),VN{x,t)  converges 
boundedly  a.e  to  functions  u(x,  t),v{x,  t),  then  the  solution  w  defined  by  w(x,  t)  =  u{x,  t)  if  — 1  <  x  <  0  and 
w{x,t)  —  v{x,  t)  if  0  <  X  <  1  converges  to  the  weak  solution  of  (1.1). 

We  will  discuss  first  the  stability  of  (5.1)-(5.2).  We  state 
Proposition  5.1.  The  boundary  operators  are  dissipative,  i.e. 


(5.5)  (ujv,P(ujv(-l,t)))jv  +  ^Ww(-l,*)^“n(-l,*)  <  0, 

(5.6)  {vN,B{vN{l,t)))j^,  -  ^v%(-l,t)AvNil,t)  <  0. 
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Proposition  5.1  implies  that  the  boundary  treatment  at  the  end-points  of  the  interval  is  stable.  Example 
of  such  operators  are  given  in  [4]. 

We  are  ready  to  state  the  stability  theorem  for  the  linear  constant  coefficient  case.  In  this  case  there 
is  no  need  for  the  spectral  viscosity  terms  and  we  will  ignore  them.  We  assume  that  /  =  Au  where  A  is 
symmetric  in  equation  (1.1)  and  in  the  same  way  /+  =  A'^u,  f~  =  A~ u  where  the  eigenvalues  of  A+  are 
nonnegative  and  those  of  A~  are  nonpositive. 

Theorem  5.2.  Let  u^,  vn  be  the  solutions  of  (5.1)-(5.2).  Define 

E{t)  =  (uA?(a:,t),Mjv(a!,t))jv  +  {«jv(a:,t),«iv(a;,t))jv, 

then 


E{t)  <  E{0) 


provided  that 
(5.7) 


n  < 


2aJo  ’ 


T2  > 


2u)o 


'^3  <  -7^^ 
2cuo 


n-Ts-  — , 
OJo 


n  > 


1 

2wo ' 
1 


T2-T4  = 


uJo 


Proof;  The  proof  follows  from  multiplying  (5.1),  (5.2)  by  and  taking  the  scalar  product.  We  use 

Proposition  5.1  and  the  following  notation 


UO  =  Un{0,  t),  Vo  =  VAr(0,  t),  Oq  =  UoA^Uo,  00  =  VqA^Vo,  Jo  = 


to  get 


(5.8) 


^  -  (n  +  r3).o-  4-  (r3  + 

1  _  _  1  _ 
+  (T2  -  ^)«0  -  (t2  +  nho  +  {t4  +  ^)/^0 


The  conditions  stated  above  for  the  n’s  guarantee  that  the  right  hand  side  of  (5.8)  is  nonpositive  and 
the  proof  is  completed. 


Remark:  The  Discontinuous  Galerkin  method  applied  to  this  problem  leads  to  the  upwinding  choice  n  = 
r4  =  0,  r2  =  T3  =  —  ^.  Another  attractive  choice  that  involves  no  splitting  of  the  fluxes  is  ti  =  r2  =  —T3  = 
~  2^- 


We  turn  now  to  the  main  purpose  of  this  section,  namely  the  proof  of  convergence  in  the  nonlinear  case 
of  the  numerical  solution  to  the  correct  entropy  solution. 

We  first  show  that  the  spectral  superviscosity  terms  do  not  create  any  problems:  consider  a  compactly 
supported  (in  [-1,1])  test  function  $(x,t)  in  (7’"[— 1,1],  m  is  to  be  specified  later. 

Lemma  5.3.  Let  (j>N^-i{x,t)  and  be  the  Legendre  interpolation  polynomials  of  in  the 

intervals  [—1,0]  (with  collocation  points  (j),  and  [0,1]  (with  collocation  points  rjj),  respectively.  Then 

(5.9)  lim  {<^;v_i,S'E(wjv))  =  0, 

N-aoo 

(5.10)  Jim  (i/;a,_i,5E(ujv))  =  0. 
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where  the  spectral  superviscosities  are  defined  in  (5.3)  and  (5.4). 


Proof:  Since  is  a  polynomial  of  degree  N  it  can  be  represented  as 

N 

Un{x)  =  ^Uk,NPki^X  +  1). 
k=0 


Therefore  from  (5.3) 


S1/(wn)  =  6 


1 

^jy2s-i 


^X{x  +  1) 

N 


dx 


S 

Un 


Y,¥{k  +  \yuk^NPk{‘ix  +  \). 

k=0 


Also  the  test  function  (j)N-i  ran  be  repre.sented  as 


N 

=  y^,^k,NPk{^X  +  1). 

Jt=0 

iFrom  the  orthogonality  of  the  Legendre  polynomials  it  follows  that 

1  ^ 

|((?i>Ar-i,SV(wAr))l  =  <^jy25-i  ^^k,NUk,Nk^ik  +  l)®(Pfc(2a;  +  l),Pfc(2x  +  1)), 

*=0 

We  can  choose  m  large  enough  such  that 

-  c 

\ik:NUk,N\  <  p- 

and  since  Uk  and  {Pki2x  +  l),P*(2a:  +  1))  are  bounded,  the  proof  is  established.  The  proof  for  (5.10)  is 
similar. 


It  is  self  evident  that  the  form  of  the  spectral  viscosity  SL  is  crucial.  In  fact  the  factor  1  -a;^  is  necessary 
in  the  proof.  Note  that 


(5.11)  {(f>N-i,SV{uN))  —  (SV{4>n-i),un)  =  {SV((Pn~i),un)n 


We  basically  proved  that  the  first  argument  in  the  scalar  product  in  the  right  hand  side  of  (5.11)  tends  to 
zero  whereas  the  second  argument  is  bounded.  The  relation  (5.11)  is  not  true  for  other  forms  of  the  spectral 
viscosity  where  the  factor  \  —  x"^  does  not  appear. 

Lemma  5.4.  Let  Tj  satisfy  (5.8),  then 


^ — )n  +  (V'A^-1,  -g^)Ar  -  {f{VN) 


(5.12) 


dt 

,  sv (ujv))  +  (^AT-I ,  SV (vjv)) 


dx 


}N 


Proof:  Taking  the  scalar  product  of  equation  (5.1)  with  4)^-1  and  (5.2)  with  ipN~i  and  denoting  by  /(O,  t)  = 
/o,  for  all  the  quantities  one  gets 


dx 


-)N  =  Ti<po 
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(5.13) 


We  integrate  now  (5.12)  with  respect  to  time  to  get 


=  -  {unH  =  =  0))  -  {vN{t  =  0),^JV-l(<  =  0)) . 


We  now  use  Lemma  3.3  in  Section  3  to  convert  wtv,  vat  to  vax,  which  are  defined  in  (1.2)  as  the 
piecewise  polynomials  having  the  values  of  un,  vn  at  the  grid  points.  Combining  then  with  Lemma  5.3  it 
follows: 

Theorem  5.5.  Let  uat  and  vn  be  the  multi-domain  Legendre  approximation  (5.1)-(5.2)  to  (1.1).  Assume 
that  the  functions  uax  and  v^x  defined  in  (1.2)  converge  boundedly  a.e.,  then  the  limit  function  is  a  weak 
solution  of  (1.1). 
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